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We present a class of solvable SO(-D) symmetric matrix models with D bosonic matrices cou- 
pled to chiral fermions. The SO(-D) symmetry is spontaneously broken due to the phase of the 
fermion integral. This demonstrates the conjectured mechanism for the dynamical generation of 
four- dimensional space-time in the TIB matrix model, which was proposed as a nonperturbative 
definition of type IIB superstring theory in ten dimensions. 



PACS immbers ll.30.Cp, ll.30.Qc, ll.25.-w 

Introduction. — One of the biggest puzzles in super- 
string theory is that the space-time dimensionality which 
naturally allows a consistent construction of the theory is 
ten instead of four. A natural resolution of this puzzle is 
to consider our 4d space time to appear dynamically and 
the other 6 dimensions to become invisible due to some 
nonperturbative effects. This may be compared to the 
situation with QCD in the early 70's, where quarks were 
introduced to explain high-energy experiments, while the 
puzzle was that none of them has ever been observed in 
reality. The understanding of quark confinement as a 
nonperturbative phenomenon in non-Abelian gauge the- 
ories was important for QCD to be recognized as the cor- 
rect theory of strong interaction. Likewise, we think it 
important to try to understand the puzzle of space-time 
dimensionality in terms of nonperturbative dynamics of 
superstring theory. Nonperturbative formulations of su- 
perstring/M theories proposed in Rcf. j^,^ may play an 
analogous role as the lattice gauge theory. 

The issue of the dynamical generation of space-time 
has been pursued ^ in the context of the IIB matrix 
model iQ] , which was proposed as a nonperturbative def- 
inition of type IIB superstring theory in ten dimensions. 
This model is a supersymmetric matrix model composed 
of 10 bosonic matrices and 16 fermionic matrices, and it 
can be thought of as the zero- volume limit of lOd SU(A^) 
super Yang-Mills theory [|lo|. The 10 bosonic matrices 
represent the dynamical space-time and the model has 
manifest SO (10) symmetry. Our four-dimensional space- 
time may be accounted for, if configurations with only 
four extended directions dominate the integration over 
the bosonic matrices. This in particular requires that 
the SO(IO) symmetry to be spontaneously broken. It is 
suggested that the phase of the fermion integral plays a 
crucial role in such a phenomenon . 

In this Letter, we present a class of solvable 80(1?) 
symmetric matrix models, where D bosonic hermitian 
matrices are interpreted as the space-time. The models 
also include chiral fermions which yield complex fermion 
integrals as in the IIB matrix model. We study the D = 4 
case explicitly and find that the SO (4) symmetry is bro- 
ken down to SO (3); namely 3 dimensional space-time is 



generated dynamically in the SO (4) invariant model. If 
we replace the fermion integral by its absolute value, the 
model is still solvable and exhibits no SSB. This result 
clearly demonstrates the conjectured mechanism for the 
dynamical generation of space-time in the IIB matrix 
model. 

The model. — The partition function of the model we 
consider is given by 
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where A^^ (/i 



1, - ■ ■ , D) are N x N hermitian matrices 
and ■0^, are TV-dimensional row and column vectors, 
respectively. (The system has an SU(Af) symmetry.) We 
assume that D is even, but we comment on a generaliza- 
tion to odd D later. The actions (||) and (||) are SO{D) 
invariant, where the bosonic matrices transform as 
a vector, and the fermion fields and -0^ transform 
as Weyl spinors. The spinor index a runs over 1, • • • ,p, 
where p = 2^/^~^ is the dimension of the spinor space. 
The pxp matrices are the gamma matrices after Weyl 
projection. The flavour index / runs over 1, ■ • ■ , Nf. We 
take the large N limit with r = Nf/N fixed (Veneziano 
limit) |15|. The fermionic part of the model can be 
thought of as the zero- volume limit of the system of Weyl 
fermions in Z3-dimensions interacting with a background 
gauge field via fundamental coupling. 

Integrating out the fermion fields, one obtains 



Z = J dyle-^'^(detl?)^^ 



(4) 



where I? is a pN x pN matrix given by T> — F^^^. In 
D — 2, we find that detl? transforms under an S0(2) 
transformation as det T> i-^ e** det V, where 6 is the angle 
of rotation [ p^. Hence the partition function (^ vanishes 
in this case |18[. In > 4, det I? is SO(-D) invariant and 
so is the model. 

The fermion determinant det T> is complex in general. 
Under parity transformation : Af — ^Ad and Af = Ai 



1 



for (i ^ D), it becomes complex conjugate. From this, 
it follows that detV becomes real if Ad = 0, or more 
generally, if n^A^ =0 for some vector n^. 

We interpret the D bosonic N x N herniitian matrices 
Af^ as the dynamical space-time as in the IIB matrix 
model 1^. The space-time has the Euclidean signature 
as a result of the Wick rotation, which is always necessary 
in path-integral formalisms. In the present model, we can 
obtain the extent of space time 



R^=( ^tr(A^)2 )^D + rp 



(5) 



using a scaling argument for arbitrary N . 

In order to probe the possible SSB of 80(1?), we first 
generalize the bosonic action as 
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5b(m) = -7V^m^tr(4) 
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We calculate the extent in the /i-th direction 
1 
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(no summation over /i) 
\nZ{m) (7) 



for arbitrary m in the large N limit. Then we take the 
limit of — > 1 (for all v) keeping the order 



TOi < TO2 < • • • < m]j 



(8) 



If do not converge to the same value, it signals the 
SSB of SO{D) symmetry 

The method. — The model (Q) with the anisotropic 
bosonic action (^ can be solved in the large limit by 
using a technique known from Random Matrix Theory 
1^ . Integrating out the bosonic matrices , 
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The normalization factor J\f in (^ is given by 

= Uim.r'/' . (12) 

Here and henceforth, we omit irrelevant m-independent 
factors in the partition function. 

The four-fcrmi action ( p^ can be written as 
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The matrix S, where we consider (a/3) and (7^) as single 
indices, is symmetric, and one can always make it real 
by choosing the representation of T ^ properly. Hence 
one can diagonalize it as 



-'a/3, 7(5 



Oal3,pT-^pTO. 



(15) 



and (|l|) can be written as 

•S'Fcrmi = ^ f ^ Oaf3,pT^ 
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aP,pT^a[3Jg 
pr Q/3 

Each square in eq. (|l^) can be linearized by a Hubbard- 
Stratonovitch transformation according to 

/2 
daexp(-£^-iQa) . (17) 

Introducing complex matrices of size Nf , we ar- 
rive at 

^ ^ " ' " ^ (18) 

(19) 

(20) 

(21) 
(22) 
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Z ^ Jf j d(TdV^dV'exp(-7VS'G -I- Sq) 

5G = Tr(at,V) ; S^^i^iM^r^ 
where the p Nf x p Nf matrix A4 is 

^ pr 

The fermionic integration yields 

where the effective action W[d-] is given by 
W[a] = 5g -IndetA^ . 



In the large N limit, the evaluation of the partition 
function amounts to solving the saddle-point equations, 
which are given by 



{apr)fg = (<rlr)fg = ^ (X -^)^^„ v^O. 
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(23) 

Assuming that the flavour SU(iV/) symmetry is not bro- 
ken, we set (Tpr = (TpT-l, where apr G C. We can further 
take apr to be real, due to ( p3| ) . Then the effective action 
reduces to 

W = Nf^^{(7prf - lndctM(cr)| , (24) 
where the p x p matrix AI (a) is given by 

Mapia) = \/2^ \/ AprOaf3,pTcrpT ■ (25) 

pr 

Thus the problem reduces to a system of finite degrees 
of freedom. 
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Exact results in 4d- — Let us solve the saddle-point 
equations explicitly in the simplest case = 4. We 
choose Fi (i = 1, 2, 3) to be Pauli matrices and r4 — il. 
The matrix Af is a 2 x 2 matrix 



M{a) 



/ Pi (Til ; 0=^^3(722 

c = ^/picri2 ; d=y/p^cr2i, 
where we have introduced the notation 



a + ib ic + d 
ic — d a — ib 

b 

d 



The saddle-point equations are 
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where A ^ a'^ + b"^ + + d"^ . Eq. (||) imphes that A 
should take one of the four possible values pi, p2, Ps 
and P4. In each case, the effective action is evaluated as 
W = Nf{l-\TiA) . 

When the parameters to obey the order (H) , the domi- 
nant saddle-point is given by A — /34. Thus the partition 
function can be obtained as 
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Using (m) we get 
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where the ± symbol should be -f for p — 1,2,3 and — 
for p = A. In the limit of rrii, — > 1 (for all v), one obtains 



Ai — A2 



A. 
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A4 



1 - r , 



which means that the SO (4) is spontaneously broken 
down to S0(3). We note that = Y^^^t^ = 4 2r 
agrees with the finite N result (|). The S_SB IS associ- 
ated with the formation of a condensate {^ai^a), which 
is invariant under S0(3), but not under full S0(4). 

The phase of the determinant. — In order to clarify 
the role played by the phase of the determinant detP, 
let us consider the model 



Z' 
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(33) 



This model can be obtained by replacing half of the Nf 
Weyl fermions ip in (0) by Weyl fermions x with opposite 
chirality. Namely, eq. (^) can be rewritten as 



Z' = I dAdi/'dV'dxdx e^(^''+^"*+^'') 



(34) 
(35) 



We use a representation of gamma matrices in which Ti 
[i — 1, • • • , (Z? — 1)) are hermitian and Yd = il- Note 
that the flavour index / now runs over / = 1, • • • , Nf/2. 

We can solve the above model with the anisotropic 
bosonic action (|^) in the large N limit using the same 
method as before. The four-fermi action reads 

'S'Fcrmi = (V'aV^^)2^a/3,75(V'7V'5) + {xlx$)^al3,'rSiXjXs) 

+ i^pixi)^c.p,ysix'^^i) + {xLrp)^cp,,s{^xi) , (36) 

where I] can be obtained from S by replacing mD hy 
—mo- Similarly to (p^), it can be diagonalized as 



-'a/3,75 



(37) 



In order to linearize (p6|), we have to introduce four sets 
of &PT- matrices, which we denote as ct^^, ct^^, (T^^ and 
(7p^ . As before, we set af^ = crf^H, where apr £ M, 



etc.. Introducing a new complex variable dp 
icrfr)^ tli6 effective action becomes 

Nf 

W' = -^{S'c-lndetM') , 



^2i< + 



M' = 



M(cr'^) M{d) 
M{a*) M{<7^) 



(38) 
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(40) 



The p X p matrices M{a^) and M{a^) are the same as 
( p5| ) except that apr is replaced by a^^ and a^^, respec- 
tively. The new p x p matrix M{a) is given by 



y.(3,pT^ pr 



(41) 



pr 



(32) The set of solutions to the saddle-point equations is richer 



than before. There are solutions with dpr — ct*^ — 0. 
In this case, the problem reduces to the previous one. 
However, there is another class of solutions in which 
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Let us consider the D 

2x2 matrix 



4 case. The matrix M is a 



M(ct) 



d + ib ic + d 
ic ~ d d — ib 

d^y/pa-ii ; 6 = VP34 0'22 

c=\/pi4 0'i2 ; d^y^p24d-21 

where we have introduced the notations 



P = ^{fnv) ^ ; Pp\ 



(42) 
(43) 
(44) 



For the first class of solutions, the effective action at 
each saddle-point is given by W = Nf{l — \np^) , where 
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V = 1,2,3,4. For the second class of solutions, the 
saddle-point equations become 



P24 0'21 



^21 = 



^12 — A Vl4 a'12 
^^■22 = A"V34 '5'22 , 



(45) 



and their complex conjugates, where A — a^+6^+c^+d^. 
Due to (^, |A| should take one of the four values p, pi4, 
P24 and P34. In each case, the effective action is evaluated 
as W^' = 7V/(l-ln|A|) . 

Thus for arbitrary m, we find that the dominant 
saddle^oint is given by the second class of the solutions 
with |A| = p and the partition function is obtained as 
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-JSlNj In p 



Using (0) we get 



1 



1 + 



(46) 



(47) 



in the limit of 1 (for all v), which means that SO (4) 

is preserved. A nonvanishing condensate (V'^Xa +XqV'£) 
breaks chiral symmetry, but not SO (4). 

Discussion. — An interesting feature of the exact re- 
sult (^) is that A4 decreases linearly as r = Nf /N is in- 
creased. At r = 1 A4 becomes zero, and the dynami- 
cal space-time becomes completely 3-dimensional. When 
r > 1, A4 becomes negative. This is possible because the 
VEV of a real positive observable is not necessarily real 
positive if the weight involves a phase. 

One can generalize the model to odd D by considering 
Dirac fermions instead of Weyl fermions. In fact, such a 
model can be obtained from the even D model considered 
here by taking the mjj — s- 00 limit. The result for the 3d 
case can thus be read off from (^) as = 1 -I- |r- for all 
fi, which preserves the SO (3) symmetry. We note that 
in the odd D models the fermion determinant for each 
flavour is real, but it is not necessarily positive. However, 
for even Nf one obtains a real positive weight, and for 
odd Nf the sign of the weight is independent of Nf . This 
explains the absence of SSB in the 3d model. 

We think that our analytical results put the conjec- 
tured mechanism for the dynamical generation of space- 
time on firmer ground. The phase of the fermion integral 
favours lower dimensional configurations and as a result 
the space-time collapses in the large N limit. On the 
other hand, the actual dimensionality of the dynamical 
space-time in the IIB matrix model is yet to be deter- 
mined. Unfortunately standard Monte Carlo simulation 
is difficult precisely due to the existence of the phase. 
However, we have recently proposed a new method to 
circumvent this problem Analytical approaches us- 
ing approximations such as the one in Ref. may also 
be useful. We hope that the models presented in this 
Letter will serve also as a testing ground for new ideas 
to understand this interesting phenomenon. 
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